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We present an update of the Fermilab-MILC Collaboration's calculation of hadronic matrix ele- 
ments for B Q -B° mixing. This work is a more extended analysis than our recent publication of the 
SU(3)-breaking ratio £ [1]. We use the asqtad staggered action for light valence quarks in com- 
bination with the Fermilab interpretation of the Sheikoleslami-Wohlert action for heavy quarks. 
The calculations use MILC's 2+1 flavor asqtad ensembles. Ensembles include four lattice spac- 
ings from approximately 0.125 fm to 0.045 fm and up/down to strange quark mass ratios as low 
as 0.05. Our calculation covers the complete set of five operators needed to describe B mixing in 
the Standard Model and beyond. In addition to an update including a fuller set of analyzed data, 
we comment on the form of the staggered XPT extrapolation function. 
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1. Introduction 

A primary focal point in the search for new physics is the precision determination of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix, which describes quark flavor physics. One process 
of great relevance in this endeavor is neutral B mixing. In the Standard Model, this mixing occurs 
via a loop process where the contributions are suppressed both by CKM matrix elements, and, for 
all but the top quark, masses of the loop quark. Therefore, new physics could present a relatively 
large signal. It has been argued that recent tension between the Standard Model and flavor physics 
experiments could be alleviated by the presence of new physics in B mixing [2, 3, 4, 5, 6, 7]. While 
the latest analyses indicate that this may not be the case [8] , the future will, no doubt, bring new 
twists, and precise calculations of the theoretical inputs to B mixing are necessary for a thorough 
understanding of quark flavor physics [9]. 

To leading order in the operator product expansion, the most general AB = 2 effective Hamil- 
tonian contains eight dimension-six operators 

5 3 

^eff=EC ; ^ + £C,4-, (1.1) 
(=1 i=l 

with 

e x = (b a Y^Lq a ) {Py^LqP), @ A = (b a Lq a ) {PRq?), 

^ 2 = {b a Lq a ){PLq^), 5 = {b a LqP) (PRq a ), (1.2) 

^3 = (b a LqP) (bl } Lq a ), 

where a and j8 are color indices, L and R are the projection operators i(l ±75), and ^1,2,3 
are obtained from ^1,2,3 by L — > R. Due to parity conservation in QCD, the matrix elements 
(B q \S'i2.3\B q ) = (Bql^i^.slBq)- This leaves five independent matrix elements. In the Standard 
Model, only matrix elements of the first three operators show up in physical observables. 
Historically, the matrix elements have been parametrized as [10, 11] 

<S°|#|*S>00 = uM 2 Bq fl q 5«Gu), (1.3) 

with coefficients c,- = ( 2 /3, ~ 5 /i2, 1/12, 1/2, 1/6). The a% are known as bag parameters and char- 
acterize the deviation from the so-called vacuum saturation-approximation [12]. Lattice QCD can 
now provide a direct calculation of the (B^|^;|S°). Because of their historical use, though, most 
lattice calculations report the bag parameters as well as the matrix elements. 

We next review the different matrix elements and combinations of matrix elements that appear 
in specific comparisons with experiment. We look at the mass difference AM q , the ratio Mi s / Mid, 
and the width difference &T q . The Standard Model expression for the B-meson mass difference 
Mi q is given by 

AM, = (^^°) m(n)\ (1.4) 

The quantities in parentheses and T}b are known factors, Vy is a CKM matrix element, and the states 
have a relativistic normalization. Measurements of Mi q have sub-percent errors [13, 14]. There- 
fore, assuming no new physics in B-mixing, constraining the CKM matrix contribution |V ; /,V,*| 
relies on our ability to determine (B q \ 6\ \B q ) precisely. 
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The ratio 



AM, 
AM,/ 



V ts 



V td 



Vu 



Mfij 
Mb, 



(1.5) 



defines Some errors cancel in this ratio, making £ a more precisely determined quantity than 
the separate matrix elements. Additionally, in CKM matrix fits, the use of t, can aid in minimizing 
correlations between lattice inputs [2]. 

The B-meson width difference can be written as [10, 11] 



Ar„ 



Gi (B°Mn)\B°\ + G 3 {B^\ff 3 (n)\B° g ) cos (j> q + 0(1 /m b ,a s 



(1.6) 



where the G;'s are comprised of known constants and short-distance coefficients, and is the CP- 
violating phase [10]. A calculation of the complete set of matrix elements allows one to have 
results not only for the two matrix elements appearing explicitly in Eq. (1.6), but also combi- 
nations of (B®\ ^1,2,3^) that are useful to phenomenologists. Specifically, AT /AM depends on 



3%6 

*>l*2>/<fl2l 



Br) and the combination R = 2 + &3 + {\/2)ff\ is useful for estimating 



1/rrib eiTors. 

Finally, including Beyond-the-Standard-Model (BSM) contributions, AM q takes the generic 
form 



am 9 = £cKm)(s°I^W 



(1.7) 



(=1 



Lattice values of the full set of matrix element, (B° q \&\\BP q ) through (B q \tff5\B q ) , are needed to 
check that a given BSM model is consistent with experiment. For examples, see Refs. [15]. 



2. The Status of Lattice Calculations 

The Fermilab Lattice and MILC Collaborations published the value of ^ = 1.268(63) in early 
2012 [1]. This joins two other published, unquenched calculations, that of the HPQCD Collab- 
oration [16] and exploratory work done by the RBC and UKQCD Collaboration [17]. Fig. 1 (a) 
compares t, from these calculations. In the Fermilab-MILC calculation, there are two dominant 
sources of error. The first is the combined error from statistics, light-quark-discretization, and 
chiral-continuum-extrapolation; this yields a 3.7% error on The second is an estimate of the ef- 
fect of wrong-spin operators that appear in the chiral extrapolation, but were not explicitly included 
in the analysis; this yields a 3.2% error on The existence of wrong-spin terms was unknown at 
the time Ref. [16] was published, and their effect was not addressed in that paper. As explained 
in Sec. 3 below, the error from "wrong-spin contributions" will not be present in our subsequent 
analyses. 

We are performing the first unquenched calculation of the full set of matrix elements. A 
quenched calculation by Becirevic et al. [18] uses static heavy valence quarks with an interpolation 
to mi,, and a linear extrapolation in the light valence mass to reach Bj. Although such approx- 
imations are no longer necessary in modern lattice-QCD calculations, only limited examples of 
improved calculations exist. For B s mixing only, a one-lattice spacing, unquenched calculation of 

^2,31-8°) has been done by the HPQCD Collaboration [19]. Several groups are actively work- 
ing on calculations of the full set of matrix elements. Our collaboration gave preliminary results 
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Figure 1: (a) ^ from three recent lattice calculations that include 2+1 sea quarks. The calculations are 
from FNAL-MILC [1], RBC-UKQCD [17], and HPQCD [16]. (b) A comparison of ensembles used in the 
calculation of Ref. [1] and the one described in these proceedings. The x and y axes show the value of 
the lattice spacing a and the sea-quark mass ratio m u / d /m s , respectively. Filled (green) circles mark the 
ensembles used in Ref. [1]. Open circles mark the ensembles used in the new analysis. The area of each 
circle is approximately proportional to the amount of data on that ensemble. The (purple) burst at the lower 
left marks the physical point. 



in the Lattice 2011 proceedings [20], and the ETM Collaboration reported on preliminary work at 
Lattice 2012 [21]. 

3. The Calculation 

We first review the differences between our published calculation of t, and our current analysis 
of the full set of matrix elements. We then discuss some details of the current analysis. 

3.1 Published versus current analysis 

Our calculation of t, in Ref [1] and the work described here and previously [20] have many 
points in common. Both use the MILC gauge configurations [22] with 2 + 1 flavors of asqtad 
staggered [23] sea quarks. To simulate the bottom valence quark, we use the Fermilab interpretation 
of the Sheikholeslami-Wohlert (clover) action [24], with the hopping parameter jq, tuned to produce 
the observed B s meson mass [25] . Light valence quarks are simulated with the asqtad action, and 
we work in the meson rest frame. 

Our on-going calculation of the full set of matrix elements improves on Ref. [1] in several 
ways. It includes two smaller lattice spacings, a pa 0.06 and 0.045 fm, for a total of four lattice 
spacings. Ensembles have sea-quark masses, with ratios of m u / d /m s = 0.1 at all lattice spacings 
and one ensemble with m u / d /m s = 0.05. On that ensemble, the light valence-quark mass range 



4 



and B^. mixing 



E.D. Freeland 



has been extended to include 0.05m s . Additionally, statistics on ensembles used in the earlier 
calculation have been increased by a factor of about four. This information is summarized in 
Fig. 1 (b). 

Our current calculation uses complete expressions for the chiral perturbation theory as de- 
scribed below and in Ref. [26]. This includes contributions from wrong-spin terms, and removes 
the need for the "wrong-spin contributions" error that appears for B, in Ref. [1]. 

In addition to the five matrix elements, our current calculation includes an evaluation of the bag 
parameters for each operator. Because our collaboration is calculating the decay constants /» [27], 
we will be able to account for correlations between the two results when extracting the these param- 
eters. Finally, combinations of the matrix elements including the ratio (B®\ ^IB®} / (B°|^i \B° q ), 
and 0r will be determined. 

3.2 Details of the new analysis 

We generate two-point meson correlation functions as well as three-point correlation functions 
for the operators in Eq. (1.2). For each meson, the two- and three-point data is fit simultaneously, 
with a common value for the B- meson energy, using constrained curve fitting [28, 29]. Simultane- 
ous fits easily account for correlations between the two- and three-point data, require fewer fits, and 
allow the two-point data to help constrain the energies, resulting in more precise fit results for the 
three -point amplitudes. Explicit operator expressions and fit functions can be found in Ref. [20]. 

With the three -point results in hand, we renormalize the matrix elements. These mix under 
renormalization, e.g. 

<^i)* = (l + a,Cii)(^i> + a*Ci2<^2). (3.1) 

The Qj are calculated using one-loop, mean-field improved lattice perturbation theory [30]. For 
a s , we use the "V" scheme [31] with four-loop running, as implemented in Ref. [32], and take 
a s = av(2/a). 

The last step in the analysis is the chiral-continuum extrapolation. We use SU(3), partially- 
quenched, heavy-meson, staggered XVY [26, 33], which generalizes the continuum calculation of 
Ref. [34] to include staggered discretization errors. With staggered light quarks, matrix elements 
of wrong-spin operators appear in the XVI [26]. These contributions vanish in the continuum, but 
can be of comparable size to other NLO chiral effects at non-zero lattice spacing. Because the five 
matrix elements (^1. 5) form a complete basis, wrong-spin contributions can be written in terms 
of them. This, fortunately, means no new low-energy constants (LECs) are introduced. Mixing 
occurs among (<^i), (^2), (^3)> an d, separately, among (^4), (^5). As an example of the mixing, 
the expression for (2?~| 0\ \B®) at next-to-leading order is 

(KMK) = h (l + + T q + Qq + f W + 

+{2p 2 + 2j3 3 )fJ b) + (2$ + 2^)0 + NLO analytic terms. (3.2) 

The terms W, T, and Q, are "correct spin" contributions from wave-function renormalization, tad- 
pole diagrams, and sunset diagrams, respectively. Terms f and Q, from tadpole and sunset dia- 
grams, are contributions from wrong-spin operators. The /3, and j8/ are the leading-order LECs for 
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the matrix element Equation (3.2) shows explicitly that no new LEC's appear in the 

chiral expression. In our analysis, we perform simultaneous fits to each group of matrix elements 
that mix. 

4. Outlook 

We have recently completed a calculation of the SU(3)-breaking ratio % used in the study of 
neutral B mixing [1]. In these proceedings, we present an update of our 2 + 1-flavor calculation of 
the five matrix elements needed to describe neutral B mixing in and beyond the Standard Model. 
This expanded calculation improves on the methods used in our calculation of t, via extended data 
sets, increased statistics, and improved analysis, particularly in the chiral-continuum extrapolation. 
Preliminary results from this analysis can be found in Ref. [20]. 
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